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The outputs of physiologic systems under neural reg- 
ulation exhibit (a) high degree of variability, (b) spacial 
and temporal fractal organization which remains invari- 
ant at different scales of observation, as well as (c) com- 
plex nonlinear properties (0, IH, |49[ ) . These inherent fea- 
tures of physiologic dynamics change significantly with 
different physiologic states such as wake and sleep, exer- 
cise and rest, circadian rhythms, as well as with patho- 
logic conditions. As different physiological states and 
pathologic perturbations correspond to changes or even 
breakdown in the mechanism of the underlying neural 
regulation, alterations in certain dynamical properties of 
physiologic signals have been found to be reliable markers 
of changes in physiologic control. 

Aging is traditionally associated with the process of 
decline of physiologic function and reduction of physio- 
logic complexity y, HO). One major hypothesis is that 
physiologic aging results from a gradual change in the 
underlying mechanisms of physiologic control — a regu- 
latory network of neural and metabolic pathways inter- 
acting through coupled cascades of nonlinear feedback 
loops on a range of time and length scales — leading 
to changes of physiologic dynamics. Under this hypothe- 
sis even ostensibly healthy elderly subjects would exhibit: 
(i) Loss of sensitivity and decreased responsiveness to ex- 
ternal and internal stimuli, leading to reduced physiologic 
variability (3); (ii) Breakdown of certain feedback loops 
acting at different time scales in the regulatory mech- 
anism of various physiologic systems. This breakdown 
would lead to loss of physiologic complexity as reflected 
in certain scale-invariant and nonlinear temporal char- 
acteristics of physiologic dynamics (^d, '47'). This hy- 
pothesis of a breakdown of physiologic complexity with 
healthy aging has recently been challenged (fTSh. Further, 
earlier studies have linked various pathologic states with 
breakdown of the scale-invariant fractal organization in 
physiologic dynamics, which is likely to result from dis- 
integration of coupled feedback loops in the regulatory 
mechanism (jH HI HI IM S US) • Thus, based on this 
hypothesis, mechanistically, physiologic processes under 
healthy aging would be categorized in the same class as 
pathologic dynamics where fractal organization and non- 
linear complexity is lost. 

A second hypothesis is that while aging may lead to re- 
duced variability, certain temporal fractal, scale- invariant 
and nonlinear structures embedded in physiologic dy- 
namics may remain unchanged. These two alternative 
hypotheses represent different notions about which as- 
pects of the physiologic control mechanisms are expected 
to change in the process of aging in contrast to the 
changes accompanying certain pathologic conditions. 

To test these two hypotheses we analyze cardiac dy- 
namics — a typical example of an output of an integrated 
physiologic system under autonomic neural regulation. 
Previous studies have shown that heart rate variability 
decreases with certain pathologic conditions (|49l . l79l). as 
well as with advanced age (STI. ItsI). Studies based on ap- 
proaches from statistical physics and nonlinear dynamics 
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revealed that heartbeat fluctuations in healthy subjects 
possess a self-similar fractal structure characterized by 
long-range power-law correlations over a range of time 
scales (|43l Isol [69l ). The scaling exponent associated 
with these power-law correlations was shown to change 
significantly with rest and exercise (fl6|. j4ll Isil ). pos- 
ture (fill. 70j_ 83). sleep and wake state (|34|)7across sleet 



stages 



Mia 



381 . 61, 62, 72) and circadian phases (j27 
|28| . 1521 ) , and to be a reliable marker of cardiac vulnera- 
bility under pathologic conditions (Q, [12, [H, [sl) . Fur- 
ther, studies have found that turbulence-like multifractal 
and nonlinear features in heartbeat dynamics are reduced 
and even lost with disease (jsol . [s^ [isl. ItH) . Several stud- 
ies have also reported reduced heart rate variability ((ttI ) 
(as also shown in Fig. [T]) , apparent loss of fractal orga- 
nization, as well as breakdown of scale-invariant corre- 
lations and certain nonlinear properties with advanced 
age (HO, [U, m, [13, [ei), suggesting that healthy aging 
is associated with changes in the neuroautonomic mech- 
anism of cardiac regulation related to disintegration of 
coupled feedback loops across a range of time scales. 

Here, we investigate how cardiac dynamics change 
with advanced age by analyzing scale-invariant, linear, 
and nonlinear characteristics of heartbeat fluctuations 
recorded from subjects during rest and sleep from two 
independent databases. 



I. DATA AND METHODS 

We analyze heartbeat interval recordings from two in- 
dependent databases. 

Fantasia Database: The Fantasia database ([63l ) con- 
tains 20 young and 20 elderly subjects. We carefully 
selected 19 healthy young subjects (9 male; 10 female) 
with an average age of 25.7 years (youngest 21; oldest 34) 
and 16 healthy elderly subjects (6 male; 10 female) with 
an average age 73.8 years (youngest 68; oldest 85). All 
subjects were recorded while watching the movie Fanta- 
sia (Disney, 1940) in a relaxed supine or semi-recumbent 
posture. These conditions were chosen to avoid the ef- 
fect which differences in the level of physical activity be- 
tween young and elderly subjects during daily routine 
might have on cardiac dynamics (Fig. [1]). The continu- 
ous ECG and respiration signals were digitized at 250 Hz. 
Each heartbeat was annotated using the ARISTOTLE 
arrhythmia detector (fssl). and each beat annotation was 
verified by visual inspection. Only intervals between two 
normal beats were considered. One young and four el- 
derly subjects (shown in Fig. [8]) were excluded from our 
analysis due to artefacts in the data. 

Sleep Heart Health Study (SHHS) Database: The Sleep 
Heart Health Study (SHHS) is a prospective cohort study 
designed to investigate the relationship between sleep 
disordered breathing and cardiovascular disease. Sub- 
jects were recorded during their habitual sleep periods of 
w 8h, and continuous ECG were recorded with 250 Hz 
(Fig. [5]) . Full details of the study design and cohort are 
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FIG. 1: Consecutive heartbeat RR intervals from a repre- 
sentative (a) young healthy and (b) elderly healthy subject 
from the Fantasia database. Under the same resting condi- 
tions elderly subjects exhibit significantly reduced heart rate 
variability. 



provided in ('461, [67l ). Details about obtaining the ECG 
and polysomnographic recordings are outlined in (68). 
Sleep apnea episodes were annotated, and heart rate data 
during apnea (obstructive and central) were excluded 
from our analysis (Fig. [2]). We selected a subset of 29 
subjects (8 males; 21 females) average age at the time 
of the first recording 75.9 years (youngest 72; oldest 84). 
The recordings were repeated 5 years later when the sub- 
jects were again screened and categorized as healthy. 
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FIG. 2: A representative subject from the SHHS database, 
(a) Consecutive heartbeat RR intervals and (b) apnea scoring. 

Detrended Fluctuation Analysis (DFA): We use the 
DFA method (fssh . which has been developed to quantify 
fractal correlations embedded in nonstationary signals, to 
estimate dynamic scale-invariant characteristics in heart- 
beat fluctuations. Compared with traditional correlation 
analyzes such as autocorrelation, power-spectrum analy- 
sis, and Hurst analysis, the advantage of the DFA method 
is that it can accurately quantify the correlation property 
of signals masked by polynomial trends, and is described 
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in detail in ([IJ, |13|, |26|, |39|, l8C 

The DFA method quantifies the detrended fluctuations 
F{n) of a signal at different time scales n. A power-law 
functional form F{n) ~ indicates the presence of self- 
similar organization in the fluctuations. The parameter 
a, called the scaling exponent, quantifies the correlation 
properties of the heartbeat signal: if a = 0.5, there is no 
correlation and the signal is white noise; if a = 1.5, the 
signal is a random walk (Brownian motion); if 0.5 < a < 
1.5, there are positive correlations, where large heartbeat 
intervals are more likely to be followed by large intervals 
(and vice versa for small heartbeat intervals). 

One advantage of the DFA method is that it can quan- 
tify signals with a > 1, which cannot be done using the 
traditional autocorrelation and R/S analyses (|l7f). In 
contrast to the conventional methods, the DFA method 
avoids spurious detection of apparent long-range correla- 
tions that are an artefact of nonstationary {73). Thus, 
the DFA method is able to detect subtle temporal struc- 
tures in highly heterogeneous physiologic time series. 

An inherent limitation of the DFA analysis is the max- 
imum time scale rimax for which the fluctuation func- 
tion F{n) can be reliably calculated. To ensure sufficient 
statistics at large scales it was shown that n,nax should 
be chosen Tirpay < N/6, where N is the length of the sig- 
nal (fisl . [26l . [80). For time scales n < rimax there is no 
bias in estimating the scaling exponent a. Thus, record- 
ings longer than 1 hour {N « 3600 beats) are sufficient 
to reliably quantify a up to time scales n — 600 beats, 
and differences in the length of the recordings between 
the Fantasia database (2 hours) and SHHS database (8 
hours) do not effect the estimate of a. Recent studies 
have tested the performance of the DFA method, when 
applied to correlated signals with patches of missing data, 
random spikes, superposed trends related to different ac- 
tivity levels and patches with different standard devia- 
tion and local correlations, as often found in heartbeat 
data dHHi). 

Both the Fantasia database and the NIH Sleep Heart 
Health Study database have used 250 Hz sampling rate 
for the EGG recordings. A precision of 0.004 sec (250 Hz) 
is more than sufficient for our analysis, since the DFA 
method as well as the MSA and FDA analyses we em- 
ploy (see below) are robust in that respect. Using a lower 
sampling rate (i.e., lower precision in the estimate of the 
RR-intervals) acts effectively as added random noise with 
an amplitude proportional to the sampling interval — in 
our case the amplitude of this "sampling noise" is more 
than two orders of magnitude smaller than the RR in- 
terval. It has been shown that adding noise with such 
a small amplitude to a fractal correlated signal does not 
effect the correlation scaling and fractal properties (flsl ). 

Magnitude and Sign Analyzes (MSA): Since the DFA 
method quantifies linear fractal characteristics related to 
two-point correlations, we have selected the MSA method 
to probe for long-term non-linear properties in the data. 
Specifically, it has been shown that signals with identical 
temporal organization, quantified by the DFA-scaling ex- 



ponent a, can exhibit very different non-linear properties 
captured by the MSA method (6). 



The MSA method (0, d) consists of the following 
steps: (i) given RRi series we obtain the increment series, 
ARRi — RRi^i — RRi] (ii) we decompose the increment 
series into a magnitude series \ARR\ and a sign series 
sign(Ai?_R); (iii) to avoid artificial trends we subtract 
the average from the magnitude series; (iv) because of 
limitations in the accuracy of the DFA method for es- 
timating the scaling exponents of anti-correlated signals 
(a < 0.5), we integrate the magnitude series (|26l): (v) we 
perform a scaling analysis using DFA; (vi) to obtain the 
scaling exponents for the magnitude series we measure 
the slope of F{n)/n on a log- log plot, where F{n) is the 
fluctuation function and n is the time scale of analysis. 



This approach is sensitive to nonlinear features in sig- 
nals ([t^ . We find that positive correlations in the mag- 
nitude series (amag > 0.5) are a reliable marker of long- 
term nonlinear properties. Thus, we employ the MSA 
as a complementary method to the DFA, because it can 
distinguish physiologic signals with identical long-range 
correlations, as quantified by the DFA method, but with 
different nonlinear properties and different temporal or- 
ganization for the sign(A_R_R) series. 



Fractal Dimension Analysis (FDA): The fractal di- 
mension D{k) is a local nonlinear measure used to quan- 
tify the irregularity of a time series (50). We estimate the 
fractal dimension using an algorithm proposed in (1231 ). 



Starting from a discrete time series, x{i), with i £ 
[1, A^], a new sparse time series is constructed in the 
following way 



x™; X (m) ,x (m + k) , . . . ,x 



N ■ 



(1) 



with TO G [1, k] where m and k are integers, and [ — 
denotes the largest integer number smaller than . 
Then a length measure for this sparse time series is de- 
fined as 



Fm{k) 




(2) 



with h = [ ^~™ J . For a time series x{i) with a fractal di- 
mension D the length Lm{k) averaged over to is a power- 
law function of the scale k: L{k) = (i(fc))„j ^ k^^. In 
the general case D can depend on the scale k. In this 
case, the local fractal dimension D{k) of the time series 
x{i) is defined as the negative local derivative of log L{k) 
as a function of log k. 
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TABLE 1: Overview of measures used. 



Abbreviation 



Measure 



Significance 



(RR) (AVNN) 
CTRR (SDNN) 

cTARR (RMSSD) 



Q,mag 



static measures 
mean of RR intervals 
std. deviation of RR 

std. deviation of ARR 

dynamic m,easures 
scaling exponent of RR 
scaling exponent of |Ai?i?| 
scaling exponent of sgn(Ai?i?) 
fractal dimension of RR 



inversely proportional to heart rate 

para- and sympathetic HRV measure sensitive 

to trends 

parasympathetic HRV measure insensitive to 
trends 

linear scale-invariant correlations 
nonlinear scale-invariant correlations 
fractal measure of directionality 
nonlinear fractal measure 



II. RESULTS 

A. Variability in heartbeat intervals and their 
increments 

We first test the possibility that advanced age in os- 
tensibly healthy subjects would lead to an increase in 
the average heart rate and to a significant reduction in 
heart rate variability — a behavior previously observed 
in subjects with congestive heart failure where under sup- 
pressed vagal tone increased heart rate is associated with 
reduced heart rate variability (figl [tqI. [Slh. We find that 
both young and elderly healthy subjects in the Fantasia 
database exhibit very similar group average interbeat in- 
tervals: {RR) ± (T = 0.9 ± 0.14 for the young group and 
{RR) ± fT = 1.06 ± 0.17 for the elderly group, where a is 
the standard deviation (Table 1^. This is in agreement 
with previous studies (.14 , .36. , 164 ). A Student's t-test 
shows no significant difference between the two groups 
with a p— value = 0.11. A very similar average heartbeat 
interval we observe for the healthy elderly subjects in the 
SHHS database with {RR) ±a ^ 0.92 ± 0.075, indicat- 
ing no significant difference (p— value = 0.07) compared 
to the group of young Fantasia subjects (Table [H). Fur- 
ther, comparing the group average heartbeat interval of 
the elderly subjects from the SHHS database with the 
same subjects recorded 5 years later, we find again no 
significant difference: {RR) ± a = 0.92 ± 0.08 at the 
first recording and {RR) ± a — 0.92 ±0.1 after 5 years 
(p— value = 0.92, Table[2]). Thus, we do not observe a sig- 
nificant change in the average heart rate with advanced 
age. 

To test whether there is a reduction in heart rate vari- 
ability with aging, we next estimate for each subject the 
standard deviation of the heartbeat intervals ctrr (of- 
ten denoted as SDNN) and the standard deviation of the 
increments in the consecutive heartbeat intervals ctarr 
(often denoted as RMSSD) (Table ^. For the young 
and elderly subjects in the Fantasia database we find a 
statistically significant difference with (i) a higher value 



for the group average (ctrr), and (ii) larger inter-subject 
variability for the young group: (o'RFi)±cr — 0.089±0.034 
for the young compared to ((Trr) ± ct = 0.051 ± 0.017 for 
the elderly subjects (p-value = 3.3 ■ 10"°*) (Table [2]). 
Similarly, we observe a significantly higher value for the 
group average (ctarr) for the young subjects in the Fan- 
tasia database ((ctarr) ± cr = 0.061 ± 0.031) compared 
to the elderly subjects (((Jarr) ± ct = 0.027 ± 0.012) 
(p— value = 9.9 • 10^^), again with a larger inter-subject 
variability for the young group (Table [5]) . We note that 
the sampling rate of 250 Hz does not effect the signifi- 
cance of the difference in (tarr between the young and el- 
derly groups, as this difference is approximately 0.034 sec, 
i.e., one magnitude larger than the sampling precission 
of 0.004 sec. 

For the group of healthy elderly subjects from the 
SHHS database we find a higher value of {o-rr) i cr = 
0.077±0.027 compared to the elderly group from the Fan- 
tasia database — a difference which could be attributed 
to the fact that the SHHS subjects were recorded during 
sleep where transitions between sleep stages are associ- 
ated with trends and larger ffuctuations in the interbeat 
interval time series (|37l , l6ll, [63). while the elderly Fanta- 
sia subjects were recorded during rest. In contrast, for 
(carr) we do not observe a significant difference between 
the elderly groups from the Fantasia and SHHS database 
(p— value = 0.74) (Table. [2]). However, we find a sig- 
nificant difference between young and elderly subjects, 
indicating a clear reduction in the heart rate variability 
with aging. 



B. Fractal Correlations 

We next test whether the temporal organization in the 
heartbeat ffuctuations changes in ostensibly healthy el- 
derly compared to young subjects. Earlier studies have 
shown that heartbeat fiuctuations exhibit self-similar 
power-law correlations over a broad range of time scales 
ranging from seconds to many hours (|43l . [69l ). and that 
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the scaling exponents associated with these power-law 
correlations change significantly with s leep and wake 
state (js^) and with pathologic conditions (|59Ll6Cll). reflect- 
ing changes in the underlying mechanism of cardiac regu- 
lation. Specifically, heartbeat fluctuations of healthy sub- 
jects during daily activity exhibit l//-like power spec- 
trum l|43l . I&ol . [69l ) with a scaling exponent awl (see Sec- 
tion During sleep this behavior changes to exponent 
a « 0.8 at time scales above 60 beats, indicating stronger 
anti-correlations in the interbeat increments ARR during 
sleep compared to wake state IH) (Fig. 3(a)). In con- 
trast, for pathologic conditions such as congestive heart 
failure earlier studies have reported a value for the expo- 
nent a closer to 1.5 — typical for random walk behavior 
(Brownian motion) and associated with loss of cardiac 
control (59 ). 

Applying the DFA method we obtain a very similar 
scaling behavior for a representative healthy young and a 
healthy elderly subject from the Fantasia database, both 
characterized by a scaling exponent a2 ~ 0.8 at interme- 



diate and large time scales (Fig. 3(b) and 3(c) ). At small 



time scales for both representative subjects we observe a 
crossover to a higher exponent of ai « 1.1 (Fig. |3(b)| 
and 3(c)). While there is certain inter-subject variabil- 



ity in the scaling functions F(n), this crossover behavior 
remains robust with a group average scaling exponent 
ai « 1.1 at small scales and a2 ~ 0.75 at large scales 
for the young subjects, and respectively ai ~ 1.2 and 
a2 « 0.8 for the elderly subjects (Appendix, Fig. [T^ . 
Our analysis indicates no significant difference in the 
scaling behavior between healthy young and healthy el- 
derly subjects under the resting conditions in the Fanta- 
sia study protocol (Table [2]). We note that our findings 
for the young and elderly Fantasia subjects (Fig. |3(b)| 
and 3(c) ) are very similar to the scaling behavior in heart- 



beat fluctuations previously reported for healthy subjects 
during sleep HJ, which exhibit a crossover from ai « 1.2 
at small time scales to a2 ~ 0.8 at intermediate and large 
time scales (Fig. |3(a) ). This similarity in the scaling 



properties of heartbeat dynamics of healthy subjects dur- 
ing sl eep (F ig. 3(a) ) and the Fantasia database subjects 
(Fig. |3(b)| and 3(c) ) may be attributed to the fact that 
under the Fantasia study protocol subjects are resting 
in a semi-recumbent/supine posture, watching a relax- 
ing movie — physiologic conditions which more closely 
resemble sleep than daytime activity. 

To confirm the validity of these findings, we further 
investigate the scale-invariant correlation properties of 
cardiac dynamics for healthy elderly subjects from the 
SHHS database, where heart rate data were recorded 
during sleep — a protocol which differs from the Fan- 
tasia study (see Section H]). In Fig. |3]wc show the DFA 
scaling curves for a representative SHHS subject with a 
crossover in the scaling behavior from ai « 1.1 at small 
time scales to a2 ~ 0.9 above 60 beats. This scaling be- 
havior is very similar to the one we find for both young 
and elderly subjects from the Fantasia database (Fig. [3]). 
Further, comparing the scaling behavior of the elderly 



subjects from the SHHS database to the same subjects 
recorded five years later, we do not find a significant dif- 
ference in the correlation scaling exponents ai and a2 
(Fig. Hand Tabled]). The results shown in Figs. [l[l]and 
in Appendix, Fig. [T2l indicate that the fractal correlation 
properties of healthy heartbeat dynamics remain stable 
and do not significantly change with advanced age. 
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FIG. 3: Fluctuation function F{n) vs. time scale n (in heart 
beat number) obtained using DFA-2 for (a) 6 h-long record 
of RR heartbeat intervals during wake and sleep from a rep- 
resentative healthy subject (MIT-BIH Normal Sinus Rhythm 
Database H^)), as well as 2 h-long records of a representa- 
tive (b) healthy young subject and (c) healthy elderly subject 
from the Fantasia database. A very similar scaling behavior 
is observed for the representative (b) young and (c) elderly 
subjects, which closely resemble the scaling behavior of the 
healthy subjects during sleep shown in (a) (MIT-BIH Nor- 
mal Sinus Rhythm Database l(63) ). indicating no change in 
the scale- invariant temporal correlations of heartbeat inter- 
vals with advanced age under healthy resting conditions. Scal- 
ing curves for all individuals are shown in Appendix, Fig. 1121 



C. Magnitude and sign scaling analysis (MSA) 

Recent studies have demonstrated that scale-invariant 
processes with identical long-range power-law correla- 
tions may be characterized by very different dynamics for 
the magnitude and sign of their fluctuations p. [37[). and 
that the information contained in the temporal organiza- 
tion of the magnitude and the sign time series is indepen- 
dent from the correlation properties of the original time 
series (i). Specifically, for cardiac dynamics of healthy 
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FIG. 4: Fluctuation function F{n) vs. time scale n (in heart 
beat number) obtained from DFA-2 for « 8 h long records 
of RR heartbeat intervals during sleep for (a) representative 
healthy elderly subject from the SHHS database and (b) the 
same elderly subject 5 years later. The very similar values for 
the exponents ai and Q2 indicate no breakdown of linear frac- 
tal correlations with advanced age under healthy conditions. 
Note the similarity with the scaling behavior for the young 
subjects, shown in Fig. [3] and Appendix, Fig. 1121 which is 
not consistent with the hypothesis of a gradual loss of scale- 
invariant complexity in the process of aging. 



subjects it was shown (jGl) that heartbeat intervals during 
routine daily activity exhibit correlation properties at in- 
termediate and large time scales characterized by scaling 
exponent a2 ~ 1, while at the same time scales the mag- 
nitude series of the increments in consecutive heartbeat 
intervals is characterized by a^^^ ~ 0.8. Further, while 
correlations reflect the linear properties of heartbeat dy- 
namics, the temporal structure of the magnitude of in- 
terbeat increments has been shown to relate to the non- 
linear properties encoded in the Fourier phases flfilzl)- 
For certain pathologic conditions such as congestive heart 
failure previous studies have reported loss of nonlinear- 
ity (|65l ) associated with a breakdown of the multifractal 
spectrum (|30l ). and reduced scaling exponent a™^^ for 
the magnitude series (Q). 

For the magnitude time series of the interbeat incre- 
ments we obtain a™**^ ~ 0-53 at small time scales and 
a™'^^ « 0.68 at intermediate and large time scales for a 
representative young subject (Fig. [Sl^a)), and very sim- 
ilar results with a^^^ w 0.53 and a™^^ w 0.72 for a 
representative elderly subject from the Fantasia database 
(Fig.inib)). The DFA scaling functions F{n) for all young 
and elderly subjects, shown in Appendix, Fig. ll3[ exhibit 
a consistent behavior among the subjects in each group 
with a smooth crossover from a group average magni- 
tude exponent a™'^^ ~ 0.53 at small and intermediate 
time scales to a'^'^^ « 0.64 at large scales for the young 
group, and a similar crossover from a group average expo- 
nent a^^^ w 0.6 at small and intermediate time scales to 
a^'^^ w 0.7 at large scales for the elderly group (Tabled]). 



To confirm these findings, we next calculate the mag- 
nitude scaling exponent of the interbeat increments for 
the elderly subjects from the SHHS database. Again we 
observe a crossover from a"^^^ w 0.52 at small scales 
to a™'^^ w 0.7 at large time scales shown in Fig. [SJc) 
for a representative elderly subject — a behavior very 
similar to the one observed for both young and elderly 
Fantasia subjects shown in Fig. [Sfa-b). Our analysis 
does not show a statistically significant difference in the 
group average magnitude scaling exponents a^^^ (with 
p-value=0.71) and a™'^^ (with p-value=0.57) between the 
elderly SHHS subjects and the elderly Fantasia subjects. 
Moreover, we find no significant difference in a'"'*^ (with 
p-value=0.24) and a™'** (with p-value=0.16) between the 
elderly SHHS subjects and the young Fantasia subjects. 

For the sign of the interbeat increments time series 
we again find no significant difference in the scaling be- 
havior between the young and elderly subjects in the 
Fantasia database with practically identical exponents 
of a^^" w 0.2 at short time scales and a^^'^ w 0.4 at 
intermediate and large time scales (Fig. [6fa-b)). A con- 
sistently similar behavior we observe for all subjects in 
the young and elderly group in the Fantasia database 
(Appendix, Fig. fH)) . where the scaling function F{n) ex- 
hibits a crossover from strongly anti-correlated behavior 
at short time scales to weaker anti-correlations at larger 
scales, respectively characterized by group average sign 
exponents a^^" w 0.24 for the young and a^^" ss 0.3 for 
the elderly subjects at small scales, and a^'^ ~ 0.47 for 
the young and a^'^ w 0.43 for the elderly subjects at 
large scales. These results indicate no significant differ- 
ence in the temporal organization of the sign series be- 
tween the young and the elderly subjects in the Fantasia 
database (Tabled. 

Repeating our sign scaling analysis for the SHHS 
database we observe a crossover from strongly anti- 
correlated behavior with an exponent a^^" w 0.2 at small 
time scales to weaker anti-correlations with a'^'^ w 0.4 at 
intermediate and large time scales, as shown in Fig. ^c- 
d). This crossover behavior is very similar to the one we 
find for both young and elderly Fantasia subjects (Fig. 
|6l[a-b) and Appendix, Fig. [T4)) . Moreover, we do not 
find a significant difference in the scaling of the sign se- 
ries for the elderly SHHS subjects and the same subjects 
five years later (Fig. EJc-d) and Tabled]). 



D. Fractal Dimension Analysis 

Finally, we employ the FDA method (see Section |l| 
to estimate the fractal dimension D{k) of a time se- 
ries ([itI [23I [soh . It has been demonstrated that the 
fractal dimension is a measure which represents the non- 
linear properties in the output of a dynamical system, 
so that two signals with identical scale-invariant corre- 
lations may be quantified by different fractal dimension 
depending on the degree of nonlinearity encoded in the 
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FIG. 5: Fluctuation function F{n) vs. time scale n (in beat 
number) obtained for the magnitude of the interbeat incre- 
ments |A_Ri?| using DFA-2 for a representative (a) healthy 
young and (b) healthy elderly subject from the Fantasia 
database, and for a representative (c) healthy elderly subject 
from the SHHS database and (d) the same subject recorded 
5 years later. All subjects exhibit a very similar scaling be- 
havior characterized by an exponent a™''^ ~ 0.7 at interme- 
diate and large time scales, very different than a™^^ = 0.5 
characteristic for linear processes with no correlations in the 
Fourier phases (05 S), which indicates that the long-term non- 
linear properties of heartbeat dynamics do not break down 
with advanced age under healthy resting conditions. This is 
in contrast to the hypothesis linking the process of healthy 
aging with a gradual loss of nonlinearity. Scaling curves for 
all individuals from the Fantasia database are shown in Ap- 
pendix, Fig. [131 



FIG. 6: Fluctuation function F(n) vs. time scale n (in 
beat numbers) obtained for the sign of the interbeat in- 
crements sign(A_R7i) using DFA-2 for a representative (a) 
healthy young and (b) healthy elderly subject in the Fan- 
tasia database, and for (c) representative healthy elderly sub- 
ject from the SHHS database and (d) the same elderly SHHS 
subject recorded 5 years later. All subjects exhibit a very 
similar scaling behavior for the sign with a crossover from 
strong anti-correlations with a^^" ~ 0.2 at small time scales 
to weaker anti-correlations with Oj^" « 0.4 at large scales, 
indicating a similar fractal organization of sympathetic and 
parasympathetic control in both young and elderly subjects 
under healthy resting conditions. Scaling curves for all in- 
dividuals in the Fantasia database are shown in Appendix, 
Fig.m 
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Fourier phases ()24l . |74| ). Our analysis shows no signifi- 
cant difference in the group average of the nonlinear frac- 
tal dimension measure D{k) between the young and the 
elderly subjects in the Fantasia database for the whole 
range of time scales except for a very short time interval 
of 6 to 8 heartbeats (Fig. 7(a) ) — time scales typical 



for sleep apnea (see Fig. |8]). At smaller and larger time 
scales the average fractal dimension D{k) converges for 
both groups (Fig. 7(a) I. Furthermore, we do not observe 



a statistically significant difference between the elderly 
subjects from the SHHS database and the same subjects 
recorded five years later (Fig. |7(b)[ ). These findings do 
not support the hypothesis that nonlinearity is reduced 
in healthy elderly subjects. 




6 8 
log,(k) 



do not find a significant difference in the scaling expo- 
nents ai and a2, characterizing the fractal scale- invariant 
temporal organization of heartbeat fluctuations, between 
young and elderly subjects (Tabled]). For the scaling 
properties of the magnitude and the sign of heartbeat 
fluctuations — which have been shown to carry addi- 
tional independent information about the nonlinear and 
linear properties of a time series (4, 6, 37[ ) — we find that 
these measures also remain unchanged when comparing 
young and healthy elderly subjects (Tabled]). Finally, for 
the fractal dimension D{k) of the heartbeat interval time 
series — an independent nonlinear measure — again con- 
trary to previous reports (f2ll ). we do not find significant 
differences between young and elderly subjects. Further- 
more, comparing longitudinal data from a group of el- 
derly subjects who were also recorded five year later, we 
find that the heart rate variability is not further reduced 
(Table 12), and that the scaling exponents qi and a2 of 
the heartbeat fluctuations, as well as the nonlinear fea- 
tures as measured by the magnitude exponent a™^^ and 
the fractal dimension D(k), remain stable. 

These findings indicate that in the process of aging the 
alterations in the underlying mechanisms of cardiac au- 
tonomic regulation are not likely to involve breakdown of 
coupling between feedback loops at different time scales 
or dominance of a particular feedback loop at a given 
time scale, as often observed with pathologic perturba- 
tions (HI, Im m. Si, m. Hi). Rather, our findings sug- 
gest a reduced refiexiveness of the neuroautonomic reg- 
ulation with aging, while the nonlinear feedback interac- 
tions across time scales between elements of the cardiac 
regulatory system remain unchanged. 



III. INTERPRETATION AND MODELING 



FIG. 7: Group average nonlinear fractal dimension D{k) 
versus time scale log2 k, where k is measured in beat num- 
bers, (a) For young and elderly healthy Fantasia subjects; (b) 
healthy elderly SHHS subjects and the same subjects 5 years 
later. There is no significant difference in the group averages 
indicated by the overlapping standard deviations except for 
the interval of scales k G [3, 6] marked in (a) by the symbol (*) 
(p— value = 1.94 x 10"'* and 6 x 10~^ correspondingly). Note 
the very similar profile of D{k) for all groups indicating no ap- 
parent loss of nonlinearity with aging, in agreement with our 
findings for the long-term nonlinear properties represented by 
the magnitude exponent a™^^ shown in Fig. [5] 



E. Summary of the results 

In agreement with previous studies (fl3 . [H, [13, [tI) we 
observe certain degree of reduction in heart rate variabil- 
ity, as measured by ctrr (SDNN) and ctarr (RMSSD), 
when comparing young to elderly subjects (Tabled]). In 
contrast to previous studies ((Sy, [13, [6^ . however, we 



Our findings indicate that scale-invariant correlation 
and nonlinear properties do not significantly change in 
healthy elderly subjects compared to young subjects. 
This is in contrast to some earlier studies, based on the 
same Fantasia database (or on a subset of it), which have 
reported loss of fractal organization in heartbeat fluctu- 
ations — a behavior resembling Brownian motion (ran- 
dom walk process) with a = 1.5 at small scales and white 
noise with a = 0.5 over large scales (|36l Ig^ ). as well as 
a significant loss of nonlinearity (^T) with healthy aging. 
A possible reason for these different findings may be the 
presence of artefacts in the data such as segments of cor- 
rupted recordings or certain periodic patterns (Fig. [5]). 
These periodic patterns strongly resemble episodes of 
sleep apnea, as shown in Fig. [9] (Top panel). Indeed, 
sleep apnea may be present in the elderly subjects from 
the Fantasia database, since they have not been specifi- 
cally screened for sleep apnea. Further, EGG recordings 
were taken when subjects were watching a calming movie 
for 2 hours in a semi-recumbent or supine posture during 
which subjects may have fallen asleep for periods of time, 
when apnea episodes are likely to occur. 
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TABLE 2: Average values and standard deviation of {RR), (trr (SDNN), ctarr (RMSSD), and DFA-2 scaling exponents 
for subjects from the Fantasia database and the SHHS database. For the Fantasia database, F{n) was fitted in the interval 
n £ [6, 16] for ai and n € [60, ^] for 02- For the SHHS database, F{n) was fitted in the interval n G [6, 16] for ai and 
n G [60,600] for 02- A two-tailed Student's t-test was performed to obtain the p-values. 
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FIG. 8: Recordings of heartbeat intervals from one young and 
four elderly subjects in the Fantasia database excluded from 
our analysis. The first four recordings contain many segments 
with well pronounced periodic patterns, the period and am- 
plitude of which are typical for sleep apnea (see Fig. [9]). The 
last recording contains segments of corrupted data. These 
artefacts strongly influence scaling, fractal and nonlinear mea- 
sures (see Figs. [Tol and llip. and can lead to spurious differ- 
ences between young and elderly subjects. 



The periodic patterns we observe in wide segments of 
the interbeat interval recordings shown in Fig. [5] and 
Fig. IHa) have a period of approximately 30 to 60 sec- 
onds, typical for apnea episodes (Fig. [9ljb)) (j3ll . [35L [66l). 
Similar apnea-like patterns are also present in the breath- 
ing records of some Fantasia subjects (Figs. [5]and[HIa)). 
These periodic patterns have a very strong effect on the 
scaling analysis, as shown in earlier studies ([26l). leading 
to a pronounced crossover at the time scale correspond- 
ing to the period of the patterns. This crossover sepa- 
rates a regime of apparent Brownian-motion-type behav- 
ior with a ~ 1.5 at smaller scales from a second regime 
of apparent white noise behavior a « 0.5 at larger scales 
(Figs. fTOFf) and fTTj) — a behavior which in earlier stud- 
ies ()3Q . |64|) has spuriously been attributed to changes in 
the cardiac neuroautonomic control due to aging. 

To model the effect which periodic patterns of sleep ap- 
nea have on the scaling properties of heartbeat intervals, 
we first generate a fractal correlated signal using the 
Makse et. al. algorithm (|48h . To account for the statisti- 
cal properties observed in heartbeat intervals, we rescale 
the signal to have the mean value {Xjj{i)) — 1, standard 
deviation ax^ = 0.05, and correlation scaling exponent 
ax^ = 0.8 (Fig. fTUfa)). which match the group mean 
(RR), standard deviation {afm) (Tabled, and scahng 
exponent value {a2) (Appendix, Fig. [12] (c, d)) of the el- 
derly subjects in the Fantasia database. To model the pe- 
riodic influence of sleep apnea on the heartbeat intervals, 
we generate a sinusoidal signal, Xs{i) = Asm {2m/T), 
with a period T = 50 (similar to the average period of 
50 heartbeats in apnea patterns) and amplitude A = 0.1 
(as observed in apnea patterns) (Fig. [TOFb)). and we su- 
perpose the sinusoidal signal Xg with the fractal cor- 
related signal Xn{i) to obtain X,-is{i) — Xj^i^i) + Xs{i) 
(Fig- Unic)). We note that X^s(i) strongly resembles the 
data shown in Fig. [8] and Fig. [9] 

Applying the DFA analysis to the fractal signal AT^ 
we obtain the scaling function -F,)(n) with a slope of 
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FIG. 9: Segments of interbeat RR interval time series for (a) 
an elderly subject from the Fantasia database excluded from 
this study (shown in Fig. [S] top panel), and (b) a subject di- 
agnosed with sleep apnea from the apnea-ECG database (|6^ ). 
Both subjects show very similar and pronounced periodic pat- 
terns with a period of about 50 beats, matching the periodic 
patterns in the breathing record (top panel (a)). These pat- 
terns strongly affect the scaling analysis as demonstrated in 
Figs. [TDl and [TTI 



0.8 across all scales — in agreement with the scaling 
exponent a — 0.8 we have found for healthy subjects 
(Fig. fTOf dj). For the sinusoidal signal Xg the scaling 
function Fs{n) exhibits a crossover at scale Ux ~ cor- 
responding to the period of Xg. For scales Ux < T, the 
fluctuation function Fs{n) exhibits an apparent scaling. 



with an exponent = 2. For scales 



71 X > T, due to the periodic property of the sinusoidal 
signal Xg, the fluctuation function Fs(n) is constant and 
independent of the scale n, i.e., Fs{n) ~ ATn"", where 
as = 0. Thus, changing the amplitude A leads to a ver- 
tical shift in Fs{n) (Fig.MLe)) IH). 

Applying the DFA analysis to our model signal X^s, 
we observe that i^,,s(n) exhibits a very pronounced kink 
(not present in i^^(n)) with a crossover at rix ~ T due to 
the sinusoidal trend (Fig. [TOT f)). The behavior of F,,s(n) 
around the kink is very similar to Fs{n) around ~ T. 
At small scales Ux < T and at large scales Ux > T the 
fluctuation function Fr^s{n) converges to the scaling be- 
havior expected for Fri{n). Testing our model for signals 
Xj^ with different values for a, we find that the position 
of the crossover tix for Fj^sip-^ does not depend on o;. 
Thus, this type of crossover behavior in the scaling for 
different subjects depends only on the period T of the 
periodic patterns embedded in the heartbeat signals. 

We find that our model in Fig. fTUTf) reproduces well 
the crossover behavior in F{n) observed for the sleep 
apnea subject (Apnea-ECG Database (fssh ) shown in 
Fig. m^b). Indeed, a very similar kink in F(n) is ob- 
served at scale n ~ 50 beats for this apnea subject, as 
shown in Fig. [TT] Moreover, we find that this behavior is 
also closely followed (as shown in Fig.fTTjl by the Fantasia 



subject in Fig. [9l[a). Adding the same sinusoidal trend 
to a real heartbeat signal from a healthy subject (MIT- 
BIH Normal Sinus Rhythm Database (|63r)) also leads to 
a very similar kink in F{rL) (Fig. Ill(). 

As we demonstrate in Fig. [Til the excluded Fantasia 
subject shown in Fig. [DJa) exhibits a scaling curve very 
similar to the curve obtained from a recording during 
sleep from a subject diagnosed with sleep apnea (Apnea- 
ECG Database ((631)). Further, our model reproduces well 
the crossover in the scaling behavior of F{n) and demon- 
strates that this crossover is due to the superposition 
of healthy heart rate dynamics and a sinusoidal trend 
with approximately the same period and amplitude as 
the periodic apnea patterns shown in Fig. [51 Our model 
reproduces also the scaling curve F{n) obtained for the 
elderly Fantasia subject excluded from this study and 
shown in Fig. [SJa), suggesting that artifacts may have 
been the reason why earlier studies (|2ll . [3^ . [6^ have re- 
ported scaling differences in heartbeat dynamics between 
young and elderly subjects. 

Our modeling results confirm that the presence of pro- 
nounced crossovers for some of the elderly subjects in the 
Fantasia database are due to periodic patterns embed- 
ded in the heart rate which strongly resemble sleep ap- 
nea episodes, and thus, cannot be attributed to changes 
in the underlying mechanism of cardiac neuroautonomic 
regulation associated with healthy aging. Since apnea is 
more prominent in elderly subjects, our modeling results 
(Figs, [ini and [TT|) explain why earlier studies using the 
same Fantasia database have reported higher values for 
the scaling exponent ai at small scales n and lower values 
for a2 at large scales n for the elderly subjects compared 
to the group of young subjects (21, 36), claiming changes 
in cardiac regulation with healthy aging. 



IV. DISCUSSION 

Our investigations demonstrate the presence of robust 
correlation, fractal and nonlinear properties in cardiac 
dynamics of healthy elderly subjects which remain sur- 
prisingly stable when compared to healthy young sub- 
jects. Specifically, we find that key dynamical character- 
istics such as the correlation scaling exponent of heart- 
beat fluctuations, the scaling exponent of the magni- 
tude and sign of interbeat increments, and the nonlin- 
ear fractal dimension measure do not signiflcantly change 
with advanced age. Because the scaling exponents a 
and the fractal dimension measure D quantify a ro- 
bust scale-invariant fractal and nonlinear structure in 
heartbeat fluctuations and have been 

shown to reflect underlying mechanisms of cardiac con- 
trol (jll, [19, 32, 33), our findings indicate that impor- 
tant aspects of heartbeat regulation do not break down 
with healthy aging. Moreover, we observe no significant 
change in these scaling and nonlinear measures when 
comparing healthy elderly subjects with the same sub- 
jects recorded five years later. 
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FIG. 10: Modeling crossover behavior in the scahng of heartbeat dynamics associated with periodic patterns, (a) Artificially 
generated fractal signal Xrj with long-range power-law correlations, average value and standard deviation as observed in healthy 
heartbeat data, (b) A sinusoidal signal Xs with period and amplitude matching the period T and amplitude A of typical sleep 
apnea patterns embedded in heartbeat interval time series as shown in Fig. [9l (c) Superposition of the signals X,^ in (a) and 
Xs in (b). Note the apparent similarity between the signal X^s and the time series shown in Fig. (9] (d) Fluctuation function 
Frfin) obtained using DFA-1 for the signal in (a), (e) Fluctuation function Fs{n) obtained using DFA-1 for the signal Xs 
in (b). The position of the crossover rix corresponds to the period T in Xs- Changing A leads to a vertical shift of Fs{n). (f) 
Fluctuation function Fr,s{n) obtained using DFA-1 for the signal Xr^s in (c). Note the appearance of a kink with a crossover at 
Tlx ~ T as observed in (e). 
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FIG. 11: Scaling functions F{n) vs. time scale n obtained 
for the heartbeat intervals using DFA-2 for a healthy subject 
taken from (MIT-BIH Normal Sinus Rhythm Database (jsl)) 
(o), a Fantasia database subject we excluded from this 
study (j63i ) (shown in Fig.[9l^a)), a subject with diagnosed sleep 
apnea (□) (shown in Fig. [9|b)), and a healthy subject with 
a superposed sinusoidal signal (V). A period of T = 50 beats 
and an amplitude of ^ = 0.1 sec were chosen for the sinu- 
soidal signal to model the effect of periodic patterns due to 
sleep apnea on the scaling function F(n). This effect leads 
to a change in the scaling exponent to a 1.5 (left of the 
crossover at T) and to q « 0.5 (right of the crossover), which 
may be the reason why earlier studies have reported loss of 
fractal organization in heartbeat fluctuations with healthy ag- 
ing QEl). 



These findings do not support the hypothesis that 
healthy aging may be associated with such a change in 
the mechanism of cardiac neuroautonomic control that 
would lead to a loss of all aspects of physiologic complex- 
ity. In contrast, we find that fundamental scale- invariant 
and nonlinear properties of heartbeat dynamics remain 
unchanged. Further, our findings do not support the hy- 
pothesis of a gradual change of cardiac dynamics under 
healthy conditions with advanced age, as key properties 
of these dynamics, including heart rate variability (Ta- 
ble [2]), remain stable in healthy elderly subjects with ad- 
vancing age . Indeed, in agreement with previous stud- 
ies (|l4l . [Sa . 113, FtsI ). we find a significant reduction in 
heart rate variability as measured by (Jrb. (SDNN) and 
o'ARR (RMSSD) (although not in the average heart rate) 
in healthy elderly subjects compared to healthy young 
subjects (Tabled]). The observed reduction in heart rate 
variability is also in agreement with decrease of the com- 
monly used Approximate Entropy (ApEn) measure with 
aging, as reported earlier 13), and often interpreted as 
loss of complexity. However, comparing elderly subjects 
with the same subjects years later we do not find a further 
reduction in interbeat variability. Moreover, we do not 
observe a loss in the scale-invariant fractal and nonlinear 
features in healthy elderly compared to healthy young 



subjects, indicating that the process of aging, even in el- 
derly healthy subjects, may not result in a gradual change 
of the mechanism of control. Our findings support the 
hypothesis that (i) only certain aspects of cardiac regu- 
lation may change with advanced age: these aspects are 
related to decreased responsiveness to external and in- 
ternal stimuli, leading to reduced heart rate variability; 
(ii) other fundamental features of the neuroautonomic 
cardiac control may remain stable and unchanged with 
healthy aging: These features are related to the network 
of nonlinear feedback loops responsible for the neuroauto- 
nomic regulation at different time scales, leadiiig to scale- 
invariant cascades in heartbeat fluctuations (|33. Issl. l45h. 

This new emerging picture of healthy aging is funda- 
mentally different from the changes in neural regulation 
of cardiac dynamics under pathologic conditions ((25l . [29l . 
[S^ . [ssh . and also differs from previous studies reporting 
breakdown of the scale-invariant and nonlinear features 
of heartbeat dynamics in elderly (|20l [2ll . [3^ |47| ) . Indeed, 
suppression of parasympathetic tone and dominance of 
sympathetic inputs, typical for subjects with congestive 
heart failure, lead to changes in cardiac dynamics associ- 
ated with higher heart rate ■ Isil). lower heart rate vari- 
ability ((tqI ). relative loss of the scale- invariant long-range 
correlations in the heartbeat fluctuations with scaling ex- 
ponent a between 1.25 and 1.4 (closer to a = 1.5 corre- 
sponding to Brownian motion, i.e., random walk) ([53), 
reduced responsiveness (Q), as well as breakdown of non- 
linearity and multifractality (fll. [sol. [3^. [65l) . In contrast 
to such pathologic perturbations, healthy aging appears 
to be accompanied only by a reduction in heart rate vari- 
ability as measured by aim and aARR, while the heart 
rate, the scaling and nonlinear properties remain on av- 
erage unchanged. This important dissociation between 
heart rate variability on one side and the scale-invariant 
and nonlinear temporal organization of heartbeat fluc- 
tuations on the other side may be specific for the pro- 
cess of aging, and suggests that the alterations in the 
cardiac control mechanism with advanced age differ con- 
ceptually from the mechanistic changes in the autonomic 
regulation associated with pathologic conditions. More 
specifically, the reduced heart rate variability with ad- 
vanced age suggests a reduced responsiveness of cardiac 
control to external and internal stimuli, and thus a re- 
duced strength of feedback interactions. However, the 
cascade of nonlinear feedback loops (fs^ . Issl . |45| ) control- 
ling the dynamics across different time scales may remain 
intact in healthy elderly subjects without breaking down 
at a particular scale or across a range of scales, as the 
scale-invariant fractal and nonlinear properties appear to 
remain stable with advanced age (Table [2]) . This is not 
the case with pathologic conditions such as congestive 
heart failure, where the self-organization of neural feed- 
back interactions indeed breaks down across time scales, 
shifting the dynamics closer to a process which is more 
random (loss of long-range power-law correlations) and 
is closer to a linear process (loss of nonlinearity and mul- 
tifractality) . 
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The value of the correlation exponent a2 « 0.8 we 
observe at intermediate and large time scales for both 
young and elderly Fantasia subjects (Figs. [3] and [4|) is 
consistent with earlier reports of a very similar value of 
a2 ~ 0.85 for healthy subjects during sleep, compared 
to a ~ 1 for the same subjects during wake and daily 
activity (j33 ). This is also in agreement with studies of 
heartbeat dynamics of healthy subjects during rest and 
exercise, with a « 0.8 for rest and a « 1.1 during ex- 
ercise (fl6l l4ll Isil ). Indeed, the Fantasia subjects were 
recorded under conditions of rest (see Data and Meth- 
ods Section IJ) ([63l). Our findings of a w 0.8 consistently 
for both healthy young and healthy elderly subjects from 
the Fantasia database are further supported by our anal- 
ysis of data from the longitudinal SHHS study, where 
the same elderly subjects were recorded during sleep sev- 
eral year later. These observations of a < 1 are not due 
to artefacts in the heartbeat time series related to sleep 
apnea, as full polysomnographic data were recorded for 
the SHHS subjects indicating the apnea episodes, and 
we have excluded the apnea segments in the data from 
our analysis. Moreover, our preliminary results (a focus 
of a subsequent study) indicate no significant differences 
between young and elderly subjects even when we ac- 
count for REM and NREM sleep stages. Since there 
is no statistically significant difference in the value of 
the scaling exponent a between the young and elderly 
subjects from both databases, the a-value lower than 1 
is not likely to be related to a mechanistic breakdown 
of cardiac control with advanced age as previously sug- 
gested (jH, HI). Rather, this decrease in a is most likely 
to be related to the normal regime of cardiac regulation 
during rest and sleep when parasympathetic tone dom- 
inates during NREM sleep stages, leading to stronger 
anti-correlations with a w 0.8 in the heartbeat fluctua- 
tions ([i3,|3l[3ll3llii|)- 

We find very similar results for the scaling expo- 
nent a™^^ for the magnitude of the interbeat incre- 
ments between young and elderly subjects in the Fan- 
tasia database (Table [2]), as well as between the young 
Fantasia subjects and the elderly subjects from the 
SHHS database (see p-values reported in Results Sec- 
tion HTC]). These findings do not support the hypotheses 
that the nonlinear properties, as measured by the mag- 
nitude scaling exponent a™'**: and encoded in the Fourier 
phases (fz^). are lost with advanced age in healthy sub- 
jects under resting conditions. We note that our re- 
sults for the magnitude exponents for the young and 
elderly subjects from both databases are in agreement 
with previous studies reporting nonlinear niagnitude cor- 
relations in healthy heartbeat dynamics ([^), and more 
specifically with the magnitude exponent values found 
in the heartbeat fluctuations of healthy subjects during 
sleep 113, 111). 

Further, as the dynamics of the sign (directionality) 
of the interbeat increments is directly related to inputs 
of the sympathetic and parasympathetic branches of the 
autonomic nervous system modulating the heart rate in 



opposite directions, our findings of similar scaling for the 
sign series for both young and elderly healthy subjects 
(Table ^ indicate that fundamental features of the car- 
diac control mechanism remain unchanged with advanced 
age. We also note that our results for the sign scaling ex- 
ponent a'^s" for the young and elderly subjects from both 
databases are in agreement with the values reported in 
previous studies for healthy subjects during rest ((4l[ ) and 
sleep (f37t ). 

While our results do not show a significant difference in 
the scaling and nonlinear properties of heartbeat dynam- 
ics between healthy young and healthy elderly subjects 
during rest and sleep, we note that under conditions of 
high levels of physical activity and stress, which are as- 
sociated with a different regime of the neuroautonomic 
control, these properties may differ between young and 
elderly subjects. 

In summary, the observations reported here do not 
support the hypothesis of a continuous gradual loss of 
the scaling and nonlinear properties of cardiac dynamics 
with advanced age under healthy conditions, as we do 
not find a statistically significant change in these prop- 
erties between the young and elderly subjects from the 
Fantasia and the SHHS databases, as well as for the el- 
derly subjects from the SHHS database and the same 
subjects recorded five years later. While cardiac dy- 
namics in healthy elderly subjects is characterized by 
markedly reduced variability compared to healthy young 
subjects, the stability we observe in key fractal and non- 
linear characteristics with advanced age does not support 
the mechanistic view of a breakdown of specific feedback 
loops at given time scales in the neuroautonomic regula- 
tion (which would lead to appearance of dominant time 
scales in the dynamics), or of a breakdown of the feed- 
back interactions in cardiac control across multiple time 
scales (which would lead to random-like behavior in the 
dynamics). Indeed, both dominant time scales and close- 
to-random behavior in cardiac dynamics, have been ob- 
served under various pathologic conditions. In contrast, 
cardiac dynamics under healthy aging appears not to be- 
long to this class of processes. Instead, our results indi- 
cate that the inherent structure and temporal organiza- 
tion in the cascades of nonlinear feedback loops under- 
lying the cardiac neuroautonomic regulation remains in- 
tact in healthy elderly subjects, thus preserving the frac- 
tal and nonlinear features in heartbeat dynamics across 
all time scales. The coupling strength of these neuronal 
feedback interactions, however, is likely to diminish with 
advanced age, leading to the observed reduction in heart 
rate variability and dampened reflexive-type responsive- 
ness in elderly compared to young healthy subjects. 
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APPENDIX A 

Results of our DFA and MSA analyses for the heart- 
beat interval recordings for all young and elderly subjects 
in the Fantasia database. 
All subjects show a consistent behavior with: 

1) A smooth crossover from ai ~ 1.1 at small time 
scales to a2 ~ 0.8 at large scales for the heartbeat in- 
tervals RR for both the young and the elderly group 

(Fig.ini). 

2) A smooth crossover from a™'*^ ~ 0.6 at small time 
scales to a™**^ « 0.7 at large scales for the magnitude of 
the interbeat increments |Ai?_R| for both the young and 
the elderly group fFig. [0)1. 

3) A crossover from a^^" « 0.3 at small time scales to 
^2^" « 0.45 at large scales for the sign of the interbeat 
increments sign(Ai?i?) for both the young and the elderly 
group fFig. [TI|) . 

The results show that these fractal correlation and 
nonlinear properties of heartbeat dynamics do not break 
down with healthy aging. 
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FIG. 12: Scaling curves F(n) versus time scale n (in beat 
numbers) obtained for the RR heartbeat intervals using 
DFA-2 for (a-b) 19 young healthy subjects and (c-d) 16 el- 
derly healthy subjects in the Fantasia database. Despite cer- 
tain inter-subject variability, there is a very common scaling 
behavior with a crossover from a higher average slope ai at 
small time scales to a lower average slope ai at large scales 
as represented by the solid lines and consistent with Fig. 
and Fig. [l] Individual curves are vertically shifted to aid vi- 
sual comparison. Group average statistics are presented in 
Table (2] Vertical dashed lines indicate the range of fit. 
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FIG. 13: Scaling curves F{n) versus time scale n (in beat 
numbers) obtained for the magnitude of the interbeat incre- 
ments \ARR\ using DFA-2 for (a-b) 19 healthy young sub- 
jects and (c-d) 16 healthy elderly subjects in the Fantasia 
database. Despite certain inter-subject variability, there is 
a common scaling behavior characterized by a group average 
exponent 02 ~ 0.7 at large scales for all groups as represented 
by the solid lines, indicating presence of long-term nonlinear 
properties encoded in the Fourier phases of the heartbeat time 
series similar to those shown in Fig. (5] Curves are vertically 
shifted for clarity. Vertical dashed lines indicate the range of 
fit. 
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FIG. 14: Scaling curves F{n) versus time scale n (in beat 
numbers) obtained for the sign time series of the interbeat in- 
crements sign(A_Ri?) using DFA-2 for (a-b) 19 healthy young 
subjects and (c-d) 16 healthy elderly subjects in the Fanta- 
sia database. All subjects exhibit a crossover from strongly 
(at small scales) to weakly (at large scales) anti-correlated 
behavior with no significant statistical difference between the 
young and elderly groups (Table [2]). Scaling curves are ver- 
tically shifted for clarity. Vertical dashed lines indicate the 
range of fit. 



